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SECTION  I 
INTRODUCTION 

The  feasibility  of  using  the  summed  output  of  a  ring  of  regu¬ 
larly-spaced  radially -oriented  horizontal  seismometers  for  estimating  the 
output  of  a  vertical  seismometer  located  at  the  center  of  the  ring  was  dis- 

j.  ^ 

cussed  in  (Potter-Roden,  1965)  for  an  assumed  noise  model  consisting  of 
isotropic  single -mode  Rayleigh  waves.  The  results  of  that  preliminary 
study  were  encouraging  enough  to  justify  a  deeper  investigation  of  the  theo¬ 
retical  noise -reduction  capabilities  of  planar  seismic  arrays  containing 
both  horizontal  and  vertical  instruments.  In  the  present  report  we  wish 
to  consider  the  general  problem  of  trying  to  use  several  concentric  rings 
of  horizontals  to  estimate  the  average  output  of  a  ring  of  verticals  (concen¬ 
tric  with  the  horizontal  rings)  in  an  isotropic  noise  field  which  may  contain 
more  than  one  mode.  The  present  problem  reduces  to  that  considered  in 
P-R  if  1)  there  is  only  one  horizontal  ring,  2)  the  vertical  ring  has  radius  0 
and  contains  only  one  seismometer,  and  3)  there  is  only  one  significant 
noise  mode. 

The  underlying  motivation  remains  the  same:  the  enhancement 
of  a  vertical  or  near-vertical  P-wave.  In  the  problem  of  trying  to  detect 
teleseismic  events,  and  also  in  much  oil  exploration  work,  the  "signal"  is 
a  vertical  P-wave  which  may  be  deeply  buried  in  surface  wave  noise.  Let 
v(t)  be  the  average  output,  due  to  surface  wave  noise,  of  a  ring  of  vertical 
seismometers,  and  let  hj^(t),  .  .  .  ,  hj^(t)  be  thf  average  outputs,  due  to 
surface  wave  noise,  of  N  rings  of  horizontal  seismometers.  The  functions 
v(t)  and  hj^(t)  will  always  be  assumed  to  be  stationary  time  series.  Since 
there  is  usually  a  degree  of  statistical  coupling  between  horizontal  and 
vertical  components  of  surface  wave  noise,  it  is  reasonable  to  try  to  design 

❖ 

Hereafter  referred  to  as  P-R 


filters  gj^{t)  to  apply  to  so  as  to  make  the  autopower  spectrum  of 


e{t)  =  v{t)  -  (gj®  hj){t)  -  ...  -  h^)  (t)  (1.  1) 

as  small  as  possible.  On  the  other  hand,  since  a  vertical  P-wave  has  no 
horizontal  component,  a  processor  of  the  type  described  by  (1.  1)  will  effect 
no  distortion  of  the  P-wave  signal. 

In  this  report  we  shall  restrict  our  attention  to  processors 
of  the  type  (1.  1).  For  any  such  processor,  we  define  the  interpolation 
error  to  be 

I.E.  =  E(f)/V(f)  (1.2) 

where 

E  =  autopower  spectrum  of  e(t) 

V  s  autopower  spectrum  of  v(t) 


If  S  s  crosspower  spectrum  between  h  and  v  and  C  .=  crosspower  spectrum 

between  h  and  h., 
n  J 


E  =  V 
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(1.3) 


where  =  Fourier  transform  of  g^.  The  interpolation  error  is  minimized  if 
the  frequency-domain  filters  G^  are  chosen  so  as  to  satisfy  the  matrix  equation 
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Filters  satisfying  (1.4)  are  called  optimum  interpolation  filters,  and  for 
such  a  set  of  filters  the  interpolation  error  (1.2)  becomes 

N 

I.E.  =  1  .  (1/v)  2  (1.5) 

nsl 

(Minimum  interpolation  error) 

All  arrays  considered  in  this  report  will  have  in  common  that 
the  seismometers  are  deployed  in  concentric  rings  of  regularly-spaced 
vertical  or  radially -oriented  horizontal  components.  The  restriction  to  this 
type  of  geometry  has  been  based  upon  several  considerations.  First,  the 
mathematical  description  of  array  response  can  be  made  more  concise  and 
understandable  if  the  radial  symmetry  of  the  assumed  noise  field  is  com¬ 
plemented  by  corresponding  symmetries  in  the  array  geometry;  the  optimum 
interpolation  filters  to  be  applied  to  the  seismometers  in  a  given  ring  are 
equal  (or  nearly  so)  for  the  arrays  and  noise  fields  which  we  shall  study, 
and  hence,  we  are  justified  in  treating  the  average  output  of  all  seismometers 
in  a  ring  as  a  single  channel.  Thus,  if  there  arc  two  rings  of  horizontals 
each  containing  six  seismometers,  we  design  only  two  filters,  not  twelve. 
Another  consideration  is  that  a  ring  average  of  regularly-spaced  vertical  or 
radially-oriented  horizontal  seismometers  is  fairly  insensitive  to  directional 
characteristics  of  the  noise  field,  provided  the  separation  between  seis¬ 
mometers  is  not  large  compared  to  wavelength;  hence,  although  we  assume 
isotropy  throughout  this  investigation,  the  answers  would  be  almost  the 
same  even  if  the  noise  field  were  assumed  to  be  nearly  unidirectional. 

The  use  of  rings  of  verticals  rather  than  a  single  central 
vertical  has  a  special  purpose:  a  vertical  ring  average  has  a  low  response 
for  certain  wavelengths.  This  property  can  be  used  to  overcome  one  of 
the  serious  problems  connected  with  horizontal-vertical  interpolation, 
namely  that  for  higher  modes  there  exist  frequencies  at  which  the  hori¬ 
zontal  motion  vanishes,  making  interpolation  of  the  vertical  component 
impossible  at  those  frequencies. 
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The  first  step  in  our  investigation,  the  derivation  of  formulas 
for  numerical  calculation  of  ring  crosspower  spectra  and  optimum  inter¬ 
polation  filters  for  a  given  array  geometry  and  multimode  surface  wave 
noise  field,  is  carried  out  in  Section  II.  These  same  quantities  may  be 
expressed  by  means  of  different  formulas,  which  are  less  useful  for  numerical 
calculation,  but  more  interesting  from  the  theoretical  point  of  view  because 
they  give  some  insight  into  the  difficult  question  of  designing  an  array  geo¬ 
metry  to  suit  a  particular  noise  field.  These  formulas  are  derived  and 
discussed  in  Section  LU.  Section  IV  describes  and  discusses  a  means  for 
graphically  representing  the  response  of  a  given  horizontal-vertical  inter¬ 
polation  array  to  surface  wave  noise  having  arbitrary  dispersion  and  hori¬ 
zontal-vertical  coupling.  Attention  is  also  directed  to  the  question  of  how 
the  array  responds  to  incident  P-waves  whose  apparent  angle  of  emergence 
is  not  quite  vertical.  Finally,  in  Section  V,  we  present  the  results  of 

numerical  calculation  of  array  response  for  several  given  array  geometries 
and  noise  fields. 
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SECTION  II 

MULTIMODE  RING  CROSSPOWER  SPECTRA  AND 
OPTIMUM  INTERPOLATION  FILTERS 


A.  CROSSPOWER  SPECTRUM  FOR  TWO  SEISMOMETERS 

Let  O  be  the  origin  of  polar  coordinate  system  in  the  plane 
of  the  earth's  surface.  The  position  of  a  point  P  is  specified  by  giving 
its  polai  coordinates  (r,  p). 

At  a  given  point  Pj,  =  (r^,  Pj^)  suppose  we  have  a  vertical 
seismometer  v^,  and  a  horizontal  seismometer  h^  which  is  radially 
oriented,  L  e. ,  the  orientation  angle  of  h^  is  p^.  At  the  point 

^2“  ^*^2*  ^2^  given  a  second  pair  of  seismometers  v  and  h  , 

where  v^  =  vertical,  h^  =  radially-oriented  horizontal. 


Now  suppose,  as  in  P-R,  section  9,  that  we  have  an 
isotropic  single -mode  Rayleigh  wave  noise  field,  where  the  energy 


propagates  across  the  plane  in  pla  .e  waves  and  is  uncorxelated  from 
direction  to  direction.  Let 


^(f)  =  inline  horizontal  autopower  spectrum 

K(f)  =  inline  horizontal -vertical  transfer  function 
k{f)  =  wavenumber  (cycles  per  unit  distance) 

We  wish  to  obtain  simple  expressions  for  V  =  crosspower 
spectrum  between  v^  and  v^,  S  =  crosspower  spectrum  between  h,  and  v  , 
and  C  =  crosspower  spectrum  between  h.  and  h.. 

^  Ct 

It  may  be  shown  (cf, ,  P-R,  section  9)  that 


C  = 


where  Uj  =  Zrtkrj,  u^  =  2TTkr2«  It  is  now  useful  to  observe  if  we  define 
F(Ui,U2)  =  j  exp^iu^cos  [G-p^]  -iu^cos  fO-pJ  jdG 

-TT 

then  formulas  (2. 1)  become 

(2.2) 

V=|k|  $-F(Uj.U2) 

a F(u., u  ) 
s  =  -iK$  •  - i— ^ 

a  Uj 

a^F(u  ,u*) 

c  =  $  .  - i — £ 

au2  auj 


,*'' - ■■■■?  'TR" 


?•  L  '  .  V  1 

,  .  V'vr  k'  ^ 

V'  .V  ^ 


Substitution  of  (2.  3)  into  the  first  formula  of  (2.  2)  yields 


V  =  |Kr$J^j(u) 


(2.4) 


In  order  to  get  the  corresponding  expressions  for  S  and  C  we  have  to 
operate  on  (2.  3)  by  taking  the  partial  derivatives  indicated  in  (2.  2). 


Using  the  standard  Bessel  function  identities 


VW=*(jn-lW  -Jn+l  W) 


(2.5) 


J  (x)  =  (-1)  J  (x) 
-n  '  '  n'  ' 


(Webster,  p.  322) 


we  find 


/“l  -U2CO»(Pi-p2)  \  _ 

=  -  ^ - 1  ■  Ji(u) 

=  -i  (uj  -UjCos  tPj  -Pj])  (jj  [u]  +  (ulj 


(2.6) 


9u  9u  _  / . 


(2.7) 
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B.  SINGLE-MODE  CROSSPOWER  SPECTRUM  FOR  TWO  RINGS  OF 
SEISMOMETERS 

We  shall  now  compute  the  crosspower  spectrum  W  for  two 

concentric  rings  of  regularly  spaced  seismometers.  Thus  let  rj^,  r^,  OCj^, 

a  be  real  numbers,  r,  s  o,  r.  ^  0,  and  let  N,  and  N  be  positive  integers. 
2  1  Z  1  & 


Sj(t)  =  (1/Nj)  ^aj(t)  +  . . .  +  and 

B^(t)  =  (l/N^)  ^b^(t)  +  ...  + 


a  (t)  is  the  output  of  a  seismometer  at  the  location  ^r^^,  —  +  (1  ^ j  and 

b  (t)  is  the  output  of  a  seismometer  at  the  location  ^r^,,  ^  ^  +  tx^j  .  It 

is  assumed  that  either  all  tlie  a  's  are  vertical  seismometers  or  they  are 

m 

all  radially -oriented  horizontals,  and  similarly  for  the  Thus  s^(t) 

is  the  average  output  of  a  ring  of  N.  regularly  spaced  seismometers  at 

distance  r.  from  the  origin,  with  the  ring  having  a  rotation  angle  (X  , 

J  ^  ^ 

j  =  1,  Z.  (Figure  Z.  1). 

We  wish  to  compute  W,  the  crosspower  spectrum  between 

s,  and  s...  If  ^  =  the  crosspower  spectrum  between  a  and  b  ,  then 

1  2  mn  m  n 


^1  ^Z 


=  -  i—  -y  y 

N,N.,  4^  ^ 


(2.9) 


1=1  n=l 


10 


ycrosspower  spec  rum  for  two  vertical  rit^g  averages) 

where  =  (u^  +  -  2ujU2C08 Uj  =  Zirkr^,  U2  =  2Trkr2. 


and  (j> 


217  m  2nn 

mn  "  Nj  “  N2  “1  “  “2 


Similarly,  if  the  a  's  are  horizontals  while  the  b  's  are 


verticals,  we  get  from  (2.  8)  that 


N,  N, 

1  iu 


W  =  S 


(2.11) 


m=l  n=l  I 

1 

•(^oPmnl  "^2  IVnjjj 


(crosspower  spectrum  for  a  horizontal  ring  average  and  a  vertical  rins 


average). 


Finally,  if  all  the  and  h  'a  are  horlzontala,  we  get 


”1  ”2 


'^  =  ‘=1.2=*'t^  Z  E  Jo  1“ 

^  ^  m=l  n=l 


(2.  12) 


-(“1 


2  X  2 

cosA  +  u_  cosf))  -2u 
1  mn  2  ^mn 


l"2)(l6  JolVn]  ■^77^2  l“r 


+  48  J4(“t 


(crosspower  spectrum  for  two  horizontal  ring  averages). 


C.  MULTIMODE  CROSSPOWER  SPECTRA 


We  shall  now  derive  the  array  crosspower  formulas  for 

multimode  noise.  We  suppose  that  the  noise  field  has  P  modes,  and 

th  * 

that  for  the  p  mode. 


p 

$  (f)  =  inline  horizontal  autopower  spectrum 
P 

K  (f)  =  inline  horizontal-vertical  transfer  function 

k  (f)  =  wavenumber  (cycles  per  unit  distance), 

for  p  =  1,2,  . .  .  ,P.  We  assume  that  any  two  modes  are  uncorrelated 

and  that  the  energy  from  each  mode  is  isotropic j  uncorrelated  from 

direction-to-direction,  and  propagates  in  plane  waves  across  the  half 

space.  Our  array  consists  of  M  rings  of  vertical  seismometers  and 

N  rings  of  radially-oriented  horizontal  seismometers.  In  the  m**^ 

vertical  ring  we  have  seismometers,  located  at  the  points  whose 

polar  coordinates  are  (d  ,  ■■  +  al.  s=12  a  .  fv,.. 

'  m’ A  m'*  *  ‘ke  average 

output  of  the  A^^  seismometers  in  the  m***  vertical  ring  is  denoted  by 
v^.  In  the  n  horizontal  ring  we  have  horizontal  seismometers, 
located  at  the  points  (r^,  ^  +  0^).  s  =  I,  2, .  .  .  B  ;  the  average  output 

of  the  B^  seismometers  in  the  n**^  horizontal  ring  is  denoted  by  h  .  The 
numbers  d^  and  are.  respectively,  the  radius  and  the  rotation  angle  ( 
them  vertical  ring,  m=  l,2,...,M;  similarly,  r  and  0  are 
respectively  the  radius  and  rotation  angle  of  the  n**^  horizontal  ring, 
n—  1|2,  N, 


In  this  report  the  symbol  "p"  appearing  as  a  superscript  is  merely  an 


index  and  is  not  an  exponent 
exponent. 


Any  other  superscripted  symbol  denotes  a 


Ihen  according  to  (2.  10)  and  our  assumption  that  there 
is  no  correlation  between  modes,  we  find  that  the  crosspower  spectrum 

is  given  by 


A 

P  1 


P  1X1 

^  n  p=l  g=i  \  J 


(2.13) 


(multimode  vertical  ring  average  crosspower  spectrum) 

where  u  P  =  2nkf(d  ^  +  d  ^  -  2d  d  coscp 

^  n  m  n  ^st 


2ns  _  2TTt 


%=  A 


a  -  a  . 
m  n 


If  m  =  n,  (2.  13)  reduces  to  a  slightly  simpler  formula  for 


the  vertical  ring  autopower  spectrum  V  =  V  : 

^  mm' 


T3  A 
P  m 


.  - 

-EE 

m  ^  8=1 


K^l  Jq  Mtt  kP.d^sin 


(2.  14) 


(multimode  vertical  ring  average  autopower  spectrum). 

Similarly,  the  crosspower  spectrum  S  between  the 
th  nm 

average  h^  of  the  n  horizontal  ring  and  the  average  v  of  the  m*^^ 

m 

vertical  ring  is  obtained  from  (2.  11); 

1  P  =a  .„p  p  I  (2-  ‘5) 

E  EE  ^ 

m  D  p=l  8=1  t=l 


.  l2Tlk^ 


^2  |u/j\ 


(multimode  crosspower  between  horizontal  and  vertical  ring  averapeal 


where  cp 


St  A 


+  a  -  9 
rn  n 


u  =  2TTk^  { r  +  d  -  2r  d  cos 
St  in  m  n  m 


Finally  we  obtain  from  (2.  12)  the  formula  for  the  crosspower 

spectrum  C  between  the  horizontal  ring  averages  h  and  h  : 

nj  ^  ®  n  j 


P  B 


1  JL  "  > 

-■nr  E  i: 

n  J  p=  1  s=  1  t^ 


.n  coscp  J 


0 


(multimode  crosspower  between  two  horizontal  ring  averages) 


u  2ns  _  2n  t 

where  =  5-  5—  +  6„  -  S. 

«  1 


U  ^  =  2nk^  (r  ^  +  r.^  -  2  r  r. 
St  \  n  J  n  j 


coscp 


(2.  16) 


(Zti  (  2r_^  r.  -  (r_^^  +  eoeip^^j 


D.  MULTICHANNEL  INTERPOLATION  FILTERS  AND  INTERPOLATION 
ERROR 

For  each  m  =1,2,...  ,M,  we  wish  to  design  a  set  of 

optimum  interpolation  filters  G,  G,.  ,  ....  G 

Im,  2m  Nm  to  be  used  in 

estimating  v  from  h,,  h_,  . . .  ,h-.  in  multimode  noise  (i.e. ,  G,  is  the 
m  i  2:  N  im 


.  r- 


■»> '.  WV'-- 
.  Y^^•  .  .  - 


*'  * 


frequency  filter  to  be  applied  to  to  eetimate  v_^,  j=l,  E, . . . ;  N). 


,::  ‘ 

'■’  isC?  0 

-  J'.. 

s  :  ■; 

'ii  '/-!’.  h 


.  •»  IH  -  7. 

According  to  {1.  4),  the  optimum  interpolation  filters 


t 


^jm  <»btained  by  solving  the  matrix  equation 


<2.17) 


Where  C_^,  and  are  given  by  (2. 15)  and  (2. 16).  It  is  interesting  to 

note  here  that  if  all  of  the  ntodal  horizontal-vertical  transfer  functions 
K  all  are  purely  imaginary  at  all  frequencies  (which  is  the  case  for 

.he  nsual  theoretical  model  of  a  layered  half-space  [Laster-Unville.  1,66]) 

tHpn  Ktr  ty  f  Ct  ^ aJ  /o  1/4  .  _  -  f 


then  by  ,2.  ,5,  ,g.  .g,, 

the  filters  Gj^  are  real  also. 


Once  the  optimum  interpolation  filters  G.  have  been 
computed,  the  optimum  interpolation  error  for  v^  is  giv”n  by  (1.5)  as 


(I.E.)  =  1 

m 


■  g  <■- 


.  S.  * 
jm  jm 


(2.  18) 
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SECTION  m 


ALTERNATIVE  CROSSPOWER  FORMULAS 

If  we  assume  a  certain  noise  field  and  specify  an  array 
geometry,  then  by  use  of  formulas  (2.  14-2.  16)  we  may  calculate  the 
array  crosspower  spectra  on  an  electronic  computer  to  any  desired 
degree  of  precision;  these  results  may  in  turn  be  used  to  calculate  the 
interpolation  filters  and  interpolation  error  (formulas  2.  17,  2. 18). 
However,  although  formulas  (2.  14  -  2.  16)  are  suitable  for  numerical 
calculation,  they  are  exceedingly  inelegant,  and  because  of  their 
complexity  they  can  not  give  insight  into  how  the  array  "works.  "  How 
many  rings  should  we  have  in  our  array?  How  does  the  response  of  a 
ring  change  if  we  double  or  triple  the  number  of  seismometers  in  the 
ring?  What  is  the  effect  of  varying  the  rotation  angles  and 
of  the  vertical  and  horizontal  rings  respectively?  For  a  theoretical,  i.  e,, 
non- computational,  study  of  questions  such  as  these,  formulas  (2.14-2.16) 
are  no  help.  It  is  of  interest,  therefore,  to  derive  alternative  expressions 
for  the  array  crosspower  spectra  which  do  give  insight  into  the  functioning 
of  the  array. 

A.  PRELIMINARIES 

The  principal  tool  in  the  derivation  of  our  alternative 
crosspower  expressions  will  be  Neumann's  formula  (Watson,  p.  358): 

/f  2  2 

U  X  +  y  .  2xy  cos  ¥ )  =  2-  e  J  (x)  J  (y)  cos 

where 

G 

q 


ri.  q  =  0 
12,  q  >  0 
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The  series  on  the  right  is  absolutely  convergent  for  all  real  jc,  y,  and 
The  convergence  is  uniform  in  x  and  y  if  x  and  y  are  restricted  to  lie  in 
any  finite  interval.  Hence  one  may  differentiate  the  series  term  by  term 
with  respect  x  or  y. 

We  shall  also  need  two  definitions  from  elementary  number 

theory.  (1)  An  integer  q  is  a  multiple  of  the  integer  N  if  and  only  if 
,  I  . 

there  Exists  an  integer  m  such  than  mN  -  q.  Thus,  for  example,  zero  is 
a  multiple  of  all  integers.  6  is  a  multiple  of  2,  7  is  not  a  multiple  of  3. 

(2)  Let  Nj  and  be  positive  integers.  The  least  common  m>iltiple  of 
Nj  and  N2  is  the  smallest  positive  integer  N  such  that,  N  is  a  multiple  of 

I  K. 

both/Nj  and  N2.  Thus,  the  least  common  multiple  of  3  and  7  is  21,  the 
least  common  multiple  of  6  and  8  is  24,  and  the  least  common  multiple 
of  5  and  10  is  10. 

Finally,  we  need  to  know  the  following  fact;  let  q  and  N 
be  integers*  with  N  ^  1,  and  let  6  be  any  real  nmtrijer.  Then 


\  ^  f‘r 


.^  =  1 


(.  Pr  ••] ) 


(3.2) 


N  cos  q  6 ,  if  q  is  a  multiple  of  N 
0  ,  if  q  is  not  a  multiple  of  N 


_  I  N  cos  q  6  ,  if  q  =  0,  ±N,  ±2N,  d:3N, 


,  if  otherwise., 


B.  INFINITE  SERIES  FORMULAS  FOR  CROSSPOWER  SPECTRA 


^  Let  Vj^  and  v^  be  averages  of  two  concentric  rings  of 

regularly  spaced  vertical  seismometers;  in  the  j  ring,  let  N.  =  number 

J 


of  seismometers,  r-  'radius,  a-  =  rotation  angle,  j  =1,2.  Then  in  a 

J  J 


aingls  mode  noise  field  (4,  K,  k),  the  crosspotrer  spectrum  Vj2  between 
Vj  and  V2  is  (2. 10) 


(3.3} 


Ki^§ 


12  ■  N,  N-  _  - 

1  2  m-1  nzrl 


A  /  ix 

2.  2  '’'o([“l  *  “2  -  ^'"1  “2  Vn]  )  ■ 

m  =  1  _ 1  '  J  A 


where 


u  =  ZTTkr,.  u,  =  ZTTkr,  ,  cp  =  ^  „ 

1  12  2  ^mn  Nj  1  2. 

Hence,  according  to  Neumann's  formula  (3.  1)  we  have 

00 

Z®n  ’^n  'In  ^  'P  * 

q  q  q  2  ^  rnn 


o  N,  00 

V  ,  =  M  *  .  V 

'  -SSI  q=0 


reversing  the  order  of  summation, 


.(3.4) 


V 


K  $ 


1,  2  N.  N 


1  2  q=0 


2J  rq'^q  -H  COS  q  cp 


m=l  n=] 


mn 


The  inner  double  summation  may  be  simplified  with  the  aid  of  (3.2). 
Summing  first  on  n,  we  see  that 


is  equal  to  0  unless  q  is  a  multiple  of  N^,  by  (3.2).  Orx  the  other  hand,  if  we 
reverse  the  order  of  summation  we  see  that  the  same  sum  is  equal  to  0 
unless  q  is  also  a  multiple  of  N^.  Thus  we  see  that  tne  inner  double  summa¬ 
tion  in  (3.  4)  is  0  unless  q  is  a  multiple  both  of  N  and  N  .  i.  e. .  unless 

1  M 

q  is  a  multiple  of  N,  where  N  =  the  least  common  multiple  of  N  and  N  . 

1  2 

So.  suppose  that  q  i8_  a  multiple  of  N.  Then,  using  (3.2)  twice,  we  get 


19 


“l  ”2 


w-n  (  fZTTir 

E  E 


m  Errm 


m=l  n=l 


‘1  -“2]) 


■S"‘  ™(fr  *“'’■•)) 

=  Nj  Nj  cos  (q  (ttj  -  a^]) 


Hence,  (3.  4)  becomes 


’^1,2 


CO 

*  E  ‘ql/qN<“l’  ^qN<“2>  >=“(-9N[<li  -Clj]  ), 


(3.6) 


\2  =|K|"*  J„(U^)J^(U3) 


Q3 

+  2  |K|  iJ2  ■^qN<“l>  ^  ■'qN<"2'  '=°“  (O^  -o^]  ) 


(single-mode  crosspower  spectrum  between  two  vertical 


ring  averages! 


where 


u.  =  Zirkr, 
J  1 


r.  =  radius, 


th 

otj  ~  rotation  angle  for  the  j  ring 

N  =  least  common  multiple  of  N  and  N 

X  2 

j  =  1,  2 

Now,  let  hj  and  be  the  averages  of  two  concentric  rings 
of  radially- oriented  horizontals,  where  u.,  r.,  a,  N.,  and  N  are  the  same’ 
as  above.  We  wish  to  derive  expressions  for  ^  crosspower  between  h 
and  and  ^^^2  “  crosspower  between  andh2.  view  of  equations  (2.2), 
(2.3),  and  (3.  1),  we  see  thatS 2  C 2  can  be  obtainedfrom  (3.  6)  asfollows 


,  -i  '•>7  '.c-  \  p 


\  q=o 


qN  qN  ‘  I'^qN  '  2'  i  ~^2' 


C  =  f  .  — 2 - 

1|  2  ^9^“, 
u  1 


‘qN  V‘“l’ V  ‘“2’““ 


t 

[qN  (aj  -  o^)]  j 


We  may  carry  out  the  indicated  differentiations  term  by  term,  to  obtain 

«  (3.7) 

Sj  2  =  iK$  Jj(Uj)  (u^)  -  2iK$  ^  J^^(U2)  cos^qN  [aj-a2j 

jsmgle-mode  crosspower  spectrum  between  horizontal  and  vertical 
ring  averages) 

(3.8) 

00 

'^1,2  °  *Ji(“i)Jj(u2)  +2J  2  J'j,  (Uj)J'j,  (a^)  c°s(flN  (aj  - 

(smgle-mode  crosspower  spectrum  between  two  horizontal  ring  averages 
We  could  rewrite  (3.  7)  and  (3.  8)  by  substituting  the  identity 

f^N<V=  ^(VlE“j]  -•’qN+l[“j]). 


Let  us  express  formulas  (3.  6-3.  8)  in  the  form 


Vj  2  =  iKp  §  •  J^(aj)J^(a^)  +  .^(N.Uj.u^) 

‘^1,2  “  ^“2^  ^  ®2  “2^ 

^1.2  "  ^“2^  ^  ®3  (N.UJ.U2) 


(3.9) 
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where 


-_ii  .■V>‘'''jiL 


ej(N,  Uj,  u^)  =  2|K 


r«E 

q= 


-v  iT . jl. 

‘  ’■  i^.^ry  'x 


^  “^qN  ''^l^  “^qN  ^*^2^  cos  JqN  (a^ -o^)],  etc. 


The  fact  that  (3.  1)  in  absolutely  convergent  implies  that 

(3. 10) 


for  any  fixed  values  of  u  and  u  we  have 

X  M 


lim  e^  (N,  Uj.u^)  =  0,  j  =  1,  2,  3. 

N-»eo 


^  ^Vf 
-  A' 


■j  ‘  ■  ii'  - 


Since  N  is  the  least  common  multiple  of  N  and  N  ,  the  numbers  of 

X  M 

seismometers  in  our  two  rings,  (3.  10)  shows  that  if  we  hold  the  radii 
of  our  rings  constant  and  let  and  N2  increase  to  infinity,  then  the 
crosspower  spectra  (3.  9)  approach  the  limiting  values 


2 

Vi,2  =  iK|  «  ■  (uj)  (u^) 

®1,2  =  ‘K»Jj(Uj)J^(a^) 

'  *•'1  S’ <“2> 

That  is,  (3. 11)  gives  the  crosspower  spectra  that  would  be  obtained  by 
using  "circular  seismometers  "  whose  outputs  are  respectively  the  average 
integrated  motion  of  all  points  in  continuous  rings  with  respective  radii  r^ 
and  r^.  It  is  unfortunate  that  such  instruments  do  not  exist,  because, 
they  would  be  the  ideal  instruments  t^'  employ  in  horizontal-vertical 
interpolation  processing. 

In  fact,  •  perfect  interpolation  would  be  possible  if  circular 
seismometers  could  be  used.  For,  suppose  our  isotropic  noise  field 


contains  P  non-zero  modes  with  modal  parameters  $  ,  K  ,  and  k 

P  P  p’ 


‘ .  ^.  ,  > '  s:  ‘>M:  :%3^tm^Mfy  'iiSi  PisI 


vertical  seismometer  v  with  radius  d,  and  P  circula;;  horizontal  ' 
seismometers  h  with  radii  r  ,  n=  1,2,  . . .  ■  P.  We  shall  show  that  h  ,  ' 

“  1  iV 

h^, . . . ,  and  h  can  be  used  to  predict  v  with  zero  interpolation  error. 

Following  our  usual  notation,  let  V  =  autopower  spectrum  of  v,  S  = 

n 

crosspower  spectrum  between  h^  and  v,  and  C  .  =  crosspower  spectrum 
between  h  and  h. .  From  (3.  1 1)  we  find 


|K  I  *  J„(2tTkd) 


Define  the  matrices  $  and  A 


§  =  diag  (5  ,  $, 


and  define  the  column  vector  K  by 


Then  the  matrix  C 
matrix  V  satisfy 


(C  ),  the  column  vector  S  =  (S  ),  and  the  l  x  1 


-,4*-c  •'";4i  -, 

•  •  •■  -i-  .-7'AV^'3.‘-<^-"'~'''..’i'^'-'--'  • 


*  *■  iS,  —  iif-‘/__T 
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Now,  the  optimum  interpolation  filters  G^ . G  are'found  by  iolving^'^i®!^ 

the  enitfifi/^n  nn  -  c  1  a\  e _ ^ 


the  equation  CG  =  S  (eq.  1.4)  for  G,  where  G  is  the  column  vector  (G  ). 
"R 11^  Htr  iri  %«f  iiA  /  Cl  1  Cl  \  1.1- 1_  A  *  m  P. 


But  by  virtue  of  (3.  13)  this  equation  becomes 


t  t 

A  «  AG  =  iA  f  K 


(3.  14) 


Assuming  the  matrix  A  to  be  nonsingular,  we  can  solve  the  equation 

A  _  1  TJ^  ^  _  t  .  • 


AG  =  iK  to  obtain 


G  =  iA"^K 


(3.  15) 


and  this  solution  is  automatically  a  solution  of  (3.  14).  (Note  that  this 
solution  18  independent  of  i,  ,  the  optimum  interpolation  filters  do 
not  depend  upon  the  autopower  spectra  of  the  various  modes  for  the  case 
of  circular  seismometers. )  The  interpolation  error  is  given  by 


1  * 

T.  E.  =  1  -  -  G  S 

j=l  •’ 


(eq.  1.5) 


(3.  16) 


Now,  regarding  1x1  matrices  as  numbers,  we  note  that 


G^Sj  *  =  (S*)‘g 


r¥  i 

_7,-. .  -  .>''.7' 


...  •  ••  .  *•..  “>  ••  ■  .  :.'r  --■  J’.'  '■  •,•■•.  •  \ 
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P.E.  =  1  -  (1/V)([?5>i«]*tl) 

=  1  -  (1/V)  (-i.fK’i'jS  A-iA“^K),  by  (3.  13  and  3.  15) 

=  1  -  (K*)^  ^  Q 

(K#)‘  $K 

Thus,  in  a  noise  field  with  P  modes,  it  is  possible  to  achieve  perfect 
estimation  of  a  circular  vertical  seismometer  by  employing  P  circular 
horizontal  seismometers. 

On  the  other  hand,  a  close  look  at  the  equations  of  the 
preceding  paragraph  will  show  that  if  we  had  tried  to  estimate  v  by  using 
than  P  circular  horizontal  seismometers,  we  should  have  obtained 
rather  poor  results  in  general.  That  is,  one  needs  at  least  P  circular 
horizontals  to  handle  P  modes  of  noise.  It  is  clear  that  this  same  rule 
of  thumb  must  apply  in  the  case  of  the  finite  arrays  which  are  the  object 
of  this  study,  since  a  circular  seismometer  is  the  ideal  limiting  case  of 
a  ring  containing  finitely  many  seismometers.  Thus,  the  most  obvious 
but  perhaps  the  most  important  design  criterion  for  horizontal-vertical 
interpolation  arrays  is 

Number  of  horizontal  rings  s  number  of  noise  modes.  (3.  17) 

We  have  seen  (eq.  3.  10)  that  as  we  let  the  number  of 
seismometers  in  each  ring  approach  infinity,  the  response  jf  an  array 
rings  of  aeismometers  approximates  the  response  of  an  array  of  circular 


These  equations  together  with  (3.  9)  show  that  the  crosspower  spectra  for  "  ■' 
finite  rings  approximate  the  crosspower  spectra  for  circular  seismometers 


i!- 


if  the  u  =  2ukr  are  small  (and  the  number  of  seismometers  in  each  hori- 

i  i  -■■■  ./■■ 


zontal  ring  is  greater  than  2).  Thus,  suppose  we  have  a  fixed  finite  ring 
array  with  each  horizontal  ring  containing  at  least  3  seismometers  and 
suppose  that  the  number  of  horizontal  rings  is  greater  than  or  equal  to  the 
number  of  noise  modes  present.  Then  the  optimum  interpolation  error 
approaches  zero  as  the  maximum  wavenumber 


ir 

'l,v 


max 
l^ps  P 


(^(«) 


approaches  zero.  Conversely  it  can  be  shown  that  the  optimum  inter¬ 
polation  error  approaches  1  as  ^  , 


min 

IspsP 


(v«) 


gets  very  large. 

We  shall  now  discuss  some  further  consequences  of  equations 
(3.  6)  “  (3.  8).  For  convenience,  let. 


(Ni,  rj,  a^)x{N^.  r^,  a^) 


♦“K 

denote  the  crosspower  between  two  rings,  where  for  the  j  ring. 


N. 

J 


=  number  of  seismometers 


-  radius 


J 


a, 

J 

j 


=  rotation  angle 


=  1,  2 
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Then  (3.  6)  -  (3.  8)  show  that 


where 


N  =  least  common  multiple  of  N,  and  N 


Thus,  for  example,  the  crosspower  spectrum  between  a  6 -element  ring  and 
^  ring  is  the  same  as  the  crosspower  spectrum  between  two 

42-element  rings.  Since  a  42-element  ring  is  a  closer  approximation  to  the 
ideal  case  of  a  circular  seismometer  than  is  either  a  6-element  ring  or  a 
7 -element  ring,  we  have  evidence  that  is  very  suggestive  of  a  possible 
advantage  to  not  having  the  same  number  of  seismometers  in  all  rings  of  tl  e 
array.  There  may  also  be  an  advantage  in  varying  the  rotation  angb^s  ct 

,  .  j 

from  rmg  to  ring.  For  example,  in  the  case  of  6-  and  7 -element  horizontal 
rings,  we  have 


In  (b)  we  have  obtained  a  much  better  approximation  to  the 
between  two  circular  seismometers  than  in  (a). 
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SECTION  IV 

ARRAY  RESPONSE  FOR  A  GIVEN  SET  OF  ’  ‘ 

INTERPOLATION  FILTERS 

A.  INTERPOLATION  ERROR  FOR  AN  ARBITRARY  NOISE  MODE 

Until  now  we  have  been  concerned  with  the  problem  of 

designing  filters  to  minimize  the  interpolation  error  for  a  given  multimode 

noise  field.  Let  us  now  take  a  slightly  different  point  of  view.  Suppose 

we  have  an  array  consisting  of  a  ring  v  of  verticals  and  N  rings  h  ,  h  , . . . , 

^N  ^tjrizontals,  and  suppose  we  bi,ve  already  designed  a  set  of  Cilters 

^2’  *  "  '  ’  is  fbe  response  of  the  given  ray  (with  the  given 

set  of  filters)  to  an  arbitrary  noise  mode  with  modal  parameters  §,  K, 

and  k?  While  the  filters  may  be  optimum  interpolation  filters  for 

some  particular  multimode  noise  field,  we  ask  how  well  the  array  will 

perform  in  noise  for  which  the  filters  may  not  be  optimum. 

For  the  noise  mode  ($,  K,k)  the  output  e(t)  of  the  Interpolation 

processor  is  given  by  (1.1).  As  before,  let  E  =  autopower  of  e,  V  *  auto-  ' 

power  of  V,  =  crosspower  between  h^  and  v,  and  C  .  =  crosspower 

between  h  and  h..  We  have  (1.  3>, 
n  1  '  ' 


N 
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E  =  V  -52  G*  S 
n  n 

-  >  G  S 
n  n 

n=l 

n=l 
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tt 

n=i  j=l 


G  C  .G. 

“  nj  j 


Now,  in  the  single-mode  crosspower  formulas  (2. 10-2.  12),  note  that  the 
only  places  where  §  and  K  occur  are  outside  the  summations  as  multi¬ 
plicative  factors. 
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We  may  therefore  write 


Ikt  ^  'aCf.k) 


S  =  iK$  -bCf.k),  and 


where  the  functions  a,  b,  and  c  are  independent  of  $  and  K.  (a,  b,  and  c 

do  depend  upon  the  array  geometry  and  upon  the  filters  G  ,  but  these  are 

n 

regarded  as  fixed.)  It  may  be  shown  that  a  and  c  are  real  and  non¬ 
negative,  and  a  ^  1;  b  may  be  complex,  but  b  is  real  if  the  filters  G  are 

n 

Substituting  (4.  1)  into  (1.3),  we  find 


Hence  the  interpolation  error  I.  E.  =  E/V'is 


V- 


Let  us  define  the  total  error  (T.  E. )  for  our  processor  to  be 


T.E.  =  E/(autopower  of  a  single  vertical  seismometer) 
=  E/|k|^§  =  a-  (I.E.) 


Then,  by  (4.  3), 


T.E.  =  a 


(4.4) 


The  total  error  T.E.  irs  a  better  measure  than  I.E.  of  how 
well  our  processor  is  succeeding  in  removing  surface-wave  noise  from  a 
vertical  P-wave  trace^,  since  T  E,  talcee  into  account  the  fact  that  some 
noise  reduction  is  effected  by  v  -venumber  aliasing  between  the  seismometers 
in  the  vertical  ring  v.  If  the  vertical  ring  has  radius  zero  and  contains  only 
one  seismometer,  then  T.  E.  =1.  E. 

B.  GRAPHICAL  REPRESENTATION  OF  ARRAY  RESPONSE 

In  the  study  of  vertical  arrays  it  has  been  found  useful  to 
represent  array  response  graphically  by  contouring  the  power  response 
of  the  array  in  the  plane  for  selected  fixed  frequencies  (see,  for 

example.  Burg,  1964,  pp.  710-712).*  A  similar  method  for  graphical 
representation  of  the  response  of  horizontal -vertical  interpolation  arr&yc 
will  be  described  in  this  subsection. 

Now,  the  response  of  a  vertical  array  to  a  wave  propagating 
across  the  surface  depends  upon  frequency,  wavenumber,  &nd  direction. 
However,  the  response  of  a  horizontal -vertical  interpolation  array  depends 
upon  yet  another  variable:  the  horizontal -vertical  transfer  function  K. 

❖  ~ 

The  term  "vertical  array"  is  used  here  to  describe  a  planar  array 
with  only  vertical -component  seismometers. 


Graphical  Representation  of  Array  Response 
for  Fixed  Frequency  f.  (Sketch) 


C.  RESPONSE  OF  ARRAY  TO  NON 


■wave 
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Suppose  that  we  have'P -waves  whose  apparent  angle  of 

,.  ■  '■  .  , .  f-.  <  .p;'':':{” 

emergence  deviates  from  the  vertical  by  an  angle  3.  As  usual,  assume 
that  the  energy  is  isotropic  with  respect  to  horizontal  azimuth  and  tm-.*  ' 
correlated  from  direction  to  direction.  Let  kgCf)  be  the  inline  wavenumber 
of  the  P -waves.  Then  there  is  an  apparent  horizontal  wavenumber,  which 


IS  approximately  given  by 


and  an  apparent  horizontal -vertical  transfer  junction  K  =  cot0.  Now  assume 
that  the  filters  are  real,  so  that  a(f,k),  b(f,k),  and  c(f,k)  are  all  real, 
where  k  is  given  by  (4.  7).  Since  K  is  also  real,  eq.  (4.4)  reduces  to 


Note  that  if  the  vertical  ring  v  has  radius  0  and  contains  only  seismometer 
then  a  =  1  and  hence  4.  8  becomes 


showing  that  in  this  case,  the  vertical  amplitude  of  incident  P -waves  is 
always  increased  by  using  the  interpolation  processor. 

For  fixed  frequencies,  it  would  be  pos<sible  to  plot  T.  E. 
as  a  jimction  of  6  and  k^,  using  (4.  8).  Although  a  detailed  investigation 
has  not  been  made,  it  appears  that  even  for  a  multi -element  vertical  ring, 
small  deviations  from  the  vertical  angle  of  emergence  will  not  significant 
ly  affect  P -wave  signal  amplitudes. 
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SECTION  V 


NUMERICAL  CALCULATIONS 


In  order  to  check  the  theory  presented  in  the  preceding 

sections  on  concrete  examples,  a  Fortran  program  for  use  on  the  IBlv 

7044  has  been  written  which  <1)  accepts  as  inputs  a)  the  guometrical 

parameters  describing  a  specific  horizontal-vertical  prediction  array 

and  b)  the  modal  parameters  K^,  for  a  specific  assumed  multi- 

mode  surface  wave  noise  field,  (2)  calculates  V  ,  S  ,  and  C  acco 

m  nm  nj 

to  (2,  14-2.  16),  (3)  calculates  the  optimum  interpolation  filters  G 

jm 

according  to  (2, 17),  and  (4)  calculates  the  optimum  interpolation  erro: 
and  total  error  according  to  (2.  20).  Provision  is  made  in  the  orocran 


I  -  the  ident  .£y  matrix.  All  computations  have  used  the  three  values 
e  =  0,  e  -  0.01,  and  e  =  0.  10  meaning,  respectively,  no  uncori'elated 
noise,  1  percent  uncorrelated  noise,  and  10  percent  uncorrelated  noise 
on  the  predicting  channels. 


In  order  to  obtain  fairly  realistic  examples  of  noise  fields, 
we  have  used  the  first  3  theoretically  computed  noise  modes  for  a  12  layer 
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CASE  A:  Modes  1  and  2.  Frequency  range:  0.  05-2  cpo 

CASE  B:  Modes  1  and  2.  F requency  range:  0.  5-1.  5  cps 

CASE  C:  Modes  1  and  2.  The  horizontal  and  vertical 
amplitudes  of  mode  2  were  multiplied  by  7, 

?st^^  to  make  the  vertical  amplitudes  of  the 
1  and  2nd  modes  approximately  equal  at  1  cps. 
Frequency  range  :  0.05-2cps. 

CASE  D:  Modes  1,2,  and  3.  Horizontal  and  vertical 
amplitudes  for  mode  2  were  again  multiplied 
by  7,  and  horizontal  and  vertical  amplitudes  for 
mode  3  were  multiplied  by  10.  Frequency  range: 
0.  5-1.  5  cps. 


«.eo,.«ca:  UBO  „o.el  ,T.  I.  Ar.ay  Kes.arch  Se„ia„„.li  Be;;!  nI/ 

IX 1-34).  The  horizontal  and  vertical  amplitude  functions  and  yP,  ”  ^ 

P  =  1,2,3  for  the  first  three  modes  are  shown  in  Figures  (5.  1  -  5.3) 

kP  in  equal  to  vP/iljP.  and  fP  is  equal  to  (uP)^  The  dispersion  curves 
for  modes  1-3  are  shown  in  Figure  (5  4). 

In  Figure  5.  2,  notice  that  the  horizontal  amplitude  for 

mode  2  vanishes  at  a  frequency  of  about  0.  5  cps.  At  this  frequency. . 

no  estimation  of  a  vertical  component  from  horizontal  components  is 

possible  for  mode  2.  Hente  by  eq.  (4.  5).  the  interpolation  error  equals  1 

for  a  noise  field  consisting  solely  trf  mode  2.  However,  by  properly 

choosing  the  radius  d  for  a  vertical  ring  v.  we  can  make  the  total  error 

small,  as  shown  by  eq.  (4.6).  For  in  eq.  (4.6).  T.E.  =  a(f .  k)  f or  H  =  O. 

and  for  the  proper  choice  of  d.  a(f.k)  can  be  made  very  small  for  the 

frequency  1.  05  cps.  The  correct  choice  for  d  turns  out  to  be  approximately 
0.  892  km. 

Four  cases  have  been  computed.  The  assumed  isotropic 
noise  fields  for  the  four  cases  were: 
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Dispersion  Curves  for  Modes  1  through  3 


The  assumed  array  geometries  for  each  case  were: 

CASE  A:  One  vertical  ring,  radius  --  0.  892  km. 

Two  horizontal  rings,  radii  =  0,  5  and  0.  25 
km.  Six  seismometers  in  each  ring, 

CASE  B:  One  vertical  ring,  radius  =  0.  0,  containing 
one  seismometer.  One  horizontal  ring, 
radius  =  0,5  km,  containing  six  seismometers. 

CASE  C:  Same  as  Case  A. 

CASE  D:  One  vertical  ring,  radius  =  0.  892  km.  Three 

horizontal  rings,  radii  =  0.5,  0.4,  and  0.25  km. 
Six  seismometers  in  each  ring. 

All  rotation  angles  were  zero  for  the  cases  considered. 

In  .Figures  5.5-5.  12,  the  computed  interpolation  filters, 
interpolation  error,  and  total  error  are  shown  for  these  four  cases,  for 
an  assumed  1  percent  uncorrelated  noise  on  the  predicting  channels. 

In  Cases  A  and  B  the  first  mode  is  much  larger  than  the 
second  mode,  so  that  we  almost  have  a  single-mode  noise  field.  However 

Figures  5.  6  and  5.  8  show  that  we  still  obtain  much  better  performance  by 
using  2  rings  of  horizontals  and  a  ring  of  verticals  than  can  be  obtained  by 
using  a  single  ring  of  norizontals  and  a  single  central  vertical. 

In  the  true  multimode  cases  C  and  D  we  have  not  obtained 
good  performance  for  frequencies  less  than  0.  8  cps,  probably  because  for 
those  frequencies  the  wavenumbers  are  small  enough  to  make  some  of  the 
horizontal  ring  crosspower  spectra  almost  as  small  as  ouar  assumed 
random  noise  level.  (By  equation  (3.  9),  the  horizontal  ring  crosspower 
spectra  approach  zero  as  k  approaches  zero.)  However,  for  frequencies 
greater  than  0.  8  cps,  useful  results  are  obtained.  It  may  ba  that  by 
varying  the  array  geometries  one  may  broaden  the  band  of  frequencies 
over  which  good  interpolation  is  possible.  Further  study  is  indicated. 


Interpolation  Errc:  and  Total  Error 
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Figure  5.7.  Computed  Interpolation  Filter  for  Case  B  {Modes  1  and  2, 
1  Percent  Uncorrelated  Noise  on  Horizontal  Ring) 
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Figure  5.8.  Interpolation  Error  (Equals  Total  Error) 
for  Case  B 


Uscor related  Noise  on  Horizontal  Rine) 


HOHHa  ^VNOIX^VHa 


FLLTER  RESPONSE 


I 

I 

1 


Figure  5.  11.  Computed  Interpolation  Filters  for  Case  D  (  1.x Mode  1, 

7  X  Mode  Z,  10  x  Mode  3,  1  Percent  Uncorrelated  Noise 
on  Horizontal  Rings) 
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Figure  5. 12.  Interpolation  Error  and  Total  Error  for  Case  D 
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CONCLUSION  :■■■■ '  ^ 

The  results  of  this  study  point  toward  important  future  i  : 

applications  of  horizontal -vertical  interpolation  arrays  in  the  fields  '  ' 

of  teleseismic  signal  detection  and  oil  exploration.  Even  in  multimode 
isotropic  Rayleigh  waves,  horizontal-vertical  interpolation  arrays 
have  the  theoretical  capability  of  canceling  almost  all  surface  wave  noise 
on  a  vertical  trace  while  passing  a  vertical  P-wave  without  distortion. 

In  Section  IV,  we  saw  that  good  performance  may  be  expected  even  if  the  ' 

P-wave  is  not  quite  vertical. 


Additional  study  is  required  of  the  problem  of  optimizing 
array  geometry.  It  appears  that  the  overall  dimensions  of  horizontal- 
vertical  prediction  arrays  will  in  general  be  smaller  than  for  vertical 
arrays  with  the  same  number  of  instruments  (all  of  the  arrays  considered 
in  Section  V  were  less  than  2  km  in  greatest  dimension).  In  fact,  were 
it  not  for  the  necessity  of  assuming  the  presence  of  uncorrelated  noise 
on  the  predicting  channels,  eqs.  (3.  18)  show  that  the  best  perfcrmance 
would  be  approached  as  the  radii  of  the  horizontal  rings  approach  zero. 
The  question  of  how  best  to  choose  the  ring  radii,  however,  has  not 
been  answered  satisfactorily.  The  same  is  true  of  the  question  of 
choosing  the  number  cf  s^iamometers  and  rotation  angles  for  each 
ring.  It  is  suggested  that  further  studies  be  undertaken  to  determine 
array  response  for  representative  noise  models,  such  as  those  in 
Section  V.  In  order  to  aid  in  the  understanding  of  the  functioning  of  a 
given  array,  it  may  be  helpful  to  prepare  array  response  plots  of  the 
type  described  in  Section  IV  B. 


Finally,  we 
one  important  shortcoming 
important  now  to  procede  f: 
to  the  experimental.  The  c 
vertical  interpolation  array 
in  the  field. 
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